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The distribution of entangled states between distant parties in an opti- 
cal network is crucial for the successful implementation of various quan- 
tum communication protocols such as quantum cryptography, telepor- 
tation and dense coding jl|, 2], S]. However, owing to the unavoidable 
loss in any real optical channel, the distribution of loss-intolerant entan- 
gled states is inevitably inflicted by decoherence, which causes a degra- 
dation of the transmitted entanglement. To combat the decoherence, 
entanglement distillation, which is the process of extracting a small set 
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arge set of less entangled states, can be 



14] . Here we report on the mesoscopic 



distillation of deterministically prepared entangled light pulses that have 
undergone non-Gaussian noise. The entangled light pulses jisl . lit] , 1?| are 
sent through a lossy channel, where the transmission is varying in time 
similarly to light propagation in the atmosphere. By employing linear op- 
tical components and global classical communication, the entanglement 
is probabilistically increased. 



Entanglement distillation has been experimentally demonstrated for spin 1/2 (or qubit) 



systems exploiting a posteriori generated polarization entangled states {(J. However, the 
implementation of a scheme that is capable of distilling entanglement of continuous variable 
systems, where information is encoded into mesoscopic carriers such as the quadratures of 
light modes, has remained an experimental challenge. It has been shown theoretically that 
if the wave function for the canonically conjugate variables of the light mode is Gaussian, 
entanglement distillation can be done only by utilizing highly non-linear (thus difficult) op- 
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(•rations no matter whet lor it is a pure or a. mixed state Jls 9]. Several protocols have boon 



put forward 10, llll, ll2|, ll3( and a proof of principle experiment on the concentration of en- 
tanglement using nonlocal and non-Gaussian operations has recently been implemented flJ ]. 

In many practical scenarios, however, the transmitted quantum state will be non- 
Gaussian: one example is the transmission of light through a turbulent atmospheric channel, 
where the attenuation coefficient will fluctuate in time, thus resulting in a non-Gaussian 
quantum state 
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191 ]. Fortunately, as we will show in this letter it is possible to distill 



entanglement that has undergone such noise by means of linear optical components, a simple 
measurement induced Gaussian operation and classical communication. 

We consider an optical field mode succinctly described by its canonically conjugated 
quadrature amplitudes, which correspond to the real and imaginary parts of the complex 
field. We denote X the amplitude quadrature and P the phase quadrature, The optical 
field can be represented by a quasi-probability distribution known as the Wigner function 
W(X, P) where X and P are eigenvalues of X and P. Having two optical fields, each 
described by the quadratures (Xa,Pa) and (Xb,Pb), the joint state can be described by 
the joint Wigner function W(Xa, Pa, Xb, Pb)- If this function is Gaussian the joint optical 
state can be fully characterized by its covariance matrix. For this case the logarithmic 
negativity (which is an entanglement monotone), denoted by LN, of the state is simply 
given by 

LN = -Log 2 fi m in, (1) 

where_jU m j„ is the smallest sympletic eigenvalue of the partial transposed covariance ma- 
trix I20J. 

Suppose now that one mode of the Gaussian entangled state is sent through a medium 
with varying attenuation. We consider iV different levels of attenuation. After the trans- 
mission the state turns into a less entangled or even unentangled state and is described by 
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FIG. 1: Schematics of the entanglement distillation protocol and the experimental 
setup, a) A weak measurement on beam B is diagnosing the state and subsequently 
used to herald the highly entangled components of the state, b) The entangled states 
are prepared by interfering two polarization squeezed beams on a 50/50 BS. The lossy 
channel is implemented by sending beam B through a controllable neutral density filter 
(ND). c) The distillation operation consists of a weak measurement (implemented 
by a 7% reflecting beam splitter) and homodyne measurement, d) The verification 
measurement employs two independent polarization detectors set to measured the 
conjugate quadratures X and P. The photocurrents are down-mixed, and digitized 
with a fast A/D converter and fed into a computer, (see method section for more 
details.) 



the convex mixture 



W(X A ,P A ,X B ,P B ) 
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PiW!(X A ,P A ,X B ,P B ) 



(2) 



where Pi is the probability for a certain transmittance and the Wigner function W( represents 
the state in channel i after transmission. The individual constituents of the mixture are all 
Gaussian functions but the sum is a non-Gaussian function. Because of this non-Gaussianity, 
distillation can be enabled solely by linear optics and feedforward as illustrated in Fig. [[J The 
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FIG. 2: (color online). Experimentally measured marginal distributions illustrating the 
effect of distillation, a) marginal distribution for the amplitude quadrature in the 
tap measurement. The vertical line indicates the threshold value chosen for this 
realization, b) Marginal distributions associated with the measurements of X and P 
of beam B (two left figures) and the joint measurements Xa + Xb and Pa — Pb (two 
right figures). The black, blue and red curves are the distribution for shot noise, the 
mixed state before distillation and after distillation, respectively. 

operation consists of a weak measurement (implemented by a 7% reflecting beam splitter 
and a homodyne detector measuring X) followed by a probabilistic heralding process where 
the remaining state is kept or discarded conditioned on the measurement outcomes. If 
the outcome of the weak measurement is larger than a specified threshold value, X t h, the 
remaining state is kept 2l|, |22J, 123, l24j, I25J, thus resulting in probabilistic recovery of the 
entanglement with a corresponding increase in LN (see method section). Note that our 
protocol cannot distill more entanglement than is contained in the most entangled component 
of the mixture in eqn. ([2]). A notable difference between our distillation approach and the 
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schemes proposed in Refs. 



26| is that our procedure relies on single copies of distributed 
entangled states whereas the protocols in Ref. 12l, |26|] are based on at least two copies. 

The experimental realization is divided into three steps: preparation, distillation, and 
verification as schematically illustrated in Fig. [TJ The entangled states are prepared by in- 
terfering two squeezed beams on a 50/50 beam splitter. The squeezed beams are generated by 
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exploiting the Kerr nonlinearity experienced by ultra-short laser pulses in optical fibers 16]. 
To ease the detection process we produce polarization squeezed states which inherently con- 
tain bright polarization components that are used as local oscillators in homodyne detection 



as described in ref. 



171 ] . Details about the generation and measurement of entanglement 



can be found in the method section. The Gaussian properties of the entangled states are 
characterized by measuring the entries of the covariance matrix, though assuming that the 
intra-correlations (such as (XaPa)) are zero due to the symmetry of the states. From the 
covariance matrix we compute the smallest sympletic eigenvalue from which we find the LN 
to be 0.76 ±0.08. 

We implement the lossy channel by inserting a neutral density filter with a variable trans- 
mittance in one of the entangled beams. The entangled beam is then transmitted through 
a channel with N = 44 different levels with corresponding transmittance from 0.1 to 1 in 
steps of 0.9/45. Combining all these realizations of the experiment a mixed state such as 
the one given by eqn. pj) is formed with the probabilities pi all being identical. However, 
after the measurement we apply an envelope function to the probabilities which allow us 
to change the probability amplitudes, thus implementing different channels. Distillation of 
entanglement is demonstrated for two different lossy channels: First we consider a discrete 
channel where the transmission randomly alternates between two different levels, and sec- 
ondly, we consider the semi-continuous channel where the transmission alternates between 
45 different levels with certain probability amplitudes. The probability distributions of the 
transmittance for the discrete channel and the continuous channel are shown in Fig. [3]-l and 
Fig. Hj-1, respectively. 

The discrete channel alternates between full transmission and 25% transmission each 
realisation occurring with a probability of 50%. After transmission in such a channel the 
resulting state is a mixture of a highly and a weakly entangled state. For this state we 
measure the Gaussian LN to be — 1.63±0.02. The Gaussian entanglement is thus completely 
lost as a result of the introduction of time-dependent loss. 

The state is then fed into the distiller and we perform homodyne measurements of beam 
A, beam B and the tap beam simultaneously. The statistics of the quadrature measurements 
on the tap beam and beam B as well as the joint distribution of Xa + Xb and Pa — Pb 
are shown in Fig. [2j From the narrowing of the joint distributions to below that of the shot 
noise, we conclude qualitatively that Gaussian entanglement has been recovered. 
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FIG. 3: (color online). Experimental and theoretical results outlining the distillation of 
an entangle state from a discrete lossy channel. The experimental results are marked 
by circles and the theoretical prediction is plotted by the red solid line. The bound 
for Gaussian entanglement is given by the blue line, and the upper bound for total 
entanglement before distillation is given by the black dashed line. Both bounds are 
surpassed by the experimental data. The weight of the two constituents in the mixed 
state after distillation for various threshold values is also experimentally investigated 
and shown in the plots labelled by 1-5. The plots explicitly show the effect of the 
distillation. The error bars depend partly on the measurement uncertainties mainly 
associated with the finite resolution of the A/D converter, and partly on the statistical 
uncertainties due to the finite measurement time and the postselection process. 

The Gaussian LN has been computed for several choices of the threshold value X t h, 
and is plotted in Fig. [3] as a function of the associated success probabilities. Furthermore, 
the probability coefficients of the two states in the mixture after distillation are shown for 
different post-selection thresholds. Note that as the threshold increases the mixture of the 
two Gaussian states reduces to a single highly entangled Gaussian state, thus demonstrating 
the act of Gaussification. Based on the experimental parameters the theoretical predictions 
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are computed and illustrated in the figure by the red curve which is seen to be in good 
agreement with the experimental results. 

The results clearly show that the amount of Gaussian entanglement is increased by the 
distillation operation. To estimate whether the total entanglement is increased, we compute 
the upper bound for the LN before distillation and verify that this bound can be surpassed 
by the Gaussian LN after distillation (see the method section). The upper bound of LN 
without the Gaussian approximation is computable from the LN of each Gaussian state in 
the mixture 20J and we find LN upper = 0.49, which is shown in Fig. [3] by the dashed black 
line. We see that for a success probability around 10~ 4 the Gaussian LN crosses the upper 
bound for entanglement, and since the state at this point is Gaussified we may conclude 
that the total entanglement of the state has indeed increased as a result of the distillation. 

Entanglement distillation also comes with a cost. As the degree of entanglement is 
increasing, the number of distilled data or, equivalently, the success probability, decreases. 
E.g. when the post selection threshold is X t h = 9 SNU, the Gaussian LN is 0.67 ± 0.09 and 
the success probability is 1.69 x 10~ 5 . The protocol has extracted only 8160 highly entangled 
states from a total set of 2.4 x 10 8 less entangled states. 

We now turn our attention to a communication channel which takes on 45 different trans- 
mission levels as opposed to the 2 level channel. The distribution of the transmittance is 
illustrated in Fig. @]-l. This channel simulates a free-space optical communications channel 
where atmospheric turbulence causes scattering and beam pointing noise [lj]]. After propa- 
gation through this channel the Gaussian LN of the mixed state is found to be —0.11 ± 0.05 
which is substantially lower than the original value of 0.76 ± 0.08. The state is subse- 
quently distilled and the change in the Gaussian LN as the threshold value increases (and 
the success probability decreases) is shown in Fig. HI We clearly see the trend that the 
entanglement available for Gaussian operations is increased, ultimately reaching the level of 
LN = 0.39 ±0.07. 

A summary of the measured values for the logarithmic negativities in the various channels 
before and after distillation is presented in table [B The demonstration of a distillation proto- 



col in this letter provides a crucial step towards the construction of a quantum repeater [27] 
for transmitting continuous variables quantum states over long distances in channels in- 
flicted by non-Gaussian noise. The various ingredients for a continuous variable quantum 
repeater that could potentially overcome non-Gaussian noise - a quantum memory 28(, a 
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FIG. 4: (color online). Experimental and theoretical results outlining the distillation of 
an entangled state from a semi-continuous lossy channel. The experimental results are 
marked by circles and the theoretical prediction by the red solid curve. The evolution 
of the mixture is directly visualized in the series of probability distributions in 1-5. 
We see that for Xth = 10 SNU the probabilities associated with low transmission levels 
are and the probability for full transmission has increased to 30% as opposed to 20% 
before distillation. It is thus clear that the highly entangled states in the mixture have 
larger weight after distillation. The theoretical prediction is given by the red curve. 



teleportation protocol [29] and an entanglement distillation protocol - have now all been 
experimentally realized and the next step is to combine some of these technologies. 



METHODS 



Theory of distillation operation 



Theoretically, the distilled state reads 



POO POO ^ 

W p (X Al P Al X B ,P B ) = / J2p i W!(X A ,P A ,X B ,P B )W (X T ,F r )dX^ T (3) 

J X th J -co i=1 
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Channels 


LN (before) 


LN (after) 


Success rate P$ 


Perfect 
Discrete 

Semi-continuous 


0.76 ±0.08 
-1.63 ±0.02 
-0.11 ±0.04 


0.76 ±0.08 
0.67 ±0.09 
0.39 ±0.07 


1 

1.69 x 10~ 5 
1.66 x 10~ 5 



TABLE I: Measured Gaussian LN before and after distillation and the corresponding 
success rates. 

where X' B = VTX B - VT^TX T , P' B = VTP B - VT^TP T , X' T = VTX T + VT^TX B 



and P' T 



TP T + y/\ — TP B , and the Wigner function W (X T , P T ) represents the vacuum 
mode entering the asymmetric beam splitter with transmittance T. The question we now 
seek to answer is whether this state is more entangled than the pre-distilled state. The above 
mentioned measure - the Gaussian LN (eqn. (OQ)) - is only valid as an entanglement monotone 
if the state is Gaussian, which is not the case for all stages of our experiment. However, 
the Gaussian LN is a good measure of entanglement useful for Gaussian operations. An 
example is continuous variable quantum teleportation where an increase in the Gaussian LN 



is directly linked with an increase in the teleportation fidelity [30] . To prove that the total 
entanglement of the state has increased as a result of distillation, we compute an upper 
bound for the LN before distillation (given by LN = Lo^l \p T \ |i where p T is the partial 
transposed density matrix of the state) and demonstrate that this bound can be surpassed 
by the Gaussian LN when the state is Gaussified after distillation (as the measure then 
becomes exact). 



Generation of entanglement 

The polarization squeezed beams are produced by launching two femto second pulses with 
balanced powers onto the two orthogonal polarization modes of two different fibers. We use 
two 13.2 m long polarization-maintaining fibers. The pump source is a Cr 4+ :YAG laser at a 
wavelength of 1500 nm, repetition rate of 163 MHz. The two orthogonal polarization modes 
in each fiber are squeezed and temporally overlapped with a relative phase of tt/2, thus 
producing a circularly polarized light beam represented by the Stokes parameter S3. The 
relative phase is achieved and controlled using an interferometric birefringence compensation 
and a locking loop based on 0.1% of the fiber output. Since the optical excitation is along 
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S 3 , the orthogonal Stokes plane spanned by the Stokes parameters Si and S 2 is "dark". 
However, it contains a vacuum squeezed state which can be easily measured in a polarization 
analyzer using the orthogonally polarized excitation. Polarization squeezing in the "dark" 
plane is thus equivalent to quadrature vacuum squeezing. We thus set S(9 sq ) — > X and 
S(9 sq + 7r/2) — > P, where S(9 sq ) and S(9 sq + n/2) are the "dark" Stokes operators of the 
squeezing and anti-squeezing directions in the dark plane. The two squeezed beams interfere 
on a 50/50 beam splitter to produce entanglement and we subsequently measure all second 
order moments between the quadratures using two homodyne detectors at 17 ± 0.5 MHz. 



Measurement of entanglement 

A general polarization, Stokes, measurement apparatus consists of a half-wave plate, a 
polarizing beam splitter (PBS) and two intensity detectors. The half-wave plate enables 
the measurement of different Stokes parameters lying in the 'dark' plane when the light 
beam is circularly (S3) polarized. The PBS outputs are measured directly by using intensity 
detectors with 98% quantum efficiency InGaAs photodiodes and with an incorporated low- 
pass filter in order to avoid ac saturation due to the laser repetition oscillation. The difference 
currents from the two detectors is produced and subsequently mixed with an electronic local 
oscillator at 17 MHz, low-pass filtered (1.9 MHz), amplified (FEMTO DHPVA-100) and 
finally digitized by an A/D exit converter (Gage CompuScope 1610) at 10 7 samples per 
second with a 16-bit resolution. After these data processing steps, the noise statistics of the 
Stokes parameters are characterized at 17 MHz relative to the optical field carrier frequency 
with a bandwidth of 1 MHz. The signal is sampled around this sideband to avoid the classical 
noise present in the frequency band around the carrier. The measurement of polarization 
entanglement is accomplished by applying identical polarization measurements on beam A 
and B with the half-wave plates set to the same angle (either measuring X or P in the 
"dark" polarization plane). For each angle, the detected photocurrent noise of beam A and 
B were simultaneously sampled 2.4 x 10 8 times and the self and cross correlations between 
the data set could easily be calculated. The covariance matrix was subsequently determined 
and the logarithmic negativity was calculated. 
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